Quantum corrections of the Dzyaloshinskii-Moriya interaction on the spin-| 
AF-Heisenberg chain in an uniform magnetic field 
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We have investigated the ground state phase diagram of the ID AF spin-i Heisenberg model 
with the staggered Dzyaloshinskii-Moriya (DM) interaction in an external uniform magnetic field 
H. We have used the exact diagonalization technique. In the absence of the uniform magnetic field 
(H = 0), we have shown that the DM interaction induces a staggered chiral phase. The staggered 
chiral phase remains stable even in the presence of the uniform magnetic field. We have identified 
that the ground state phase diagram consists of four Luttinger liquid, staggered chiral, spin-flop, 
and ferromagnetic phases. 
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I. INTRODUCTION 



The effect of an external magnetic field on the quan- 
tum properties of the ID antiferromagnetic (AF) spin-i 
model has attracted much interest in recent years. Ex- 
perimental and theoretical studies of this system have 
revealed a plethora of quantum fluctuation phenomena, 
not usually observed in higher dimensions. The Hamil- 
tonian of this model in a uniform magnetic field (H) on 
a periodic chain of N sites is given by 



H - y^\JSj.~§ j+ i - HSJ], 



(1) 



where J > is the exchange coupling and H is a uniform 
magnetic field. Theoretically, in the absence of the exter- 
nal magnetic field, H = 0, the exact solution is given by 
the Bethe ansatai. The spectrum is gapless and in the 
ground state, the system is in the Luttinger liquid phase, 
where the decay of correlations follow a power law. When 
a uniform magnetic field is applied the spectrum of the 
system remains gapless until the critical field H c = 2 J. 
Here a phase transition of the Pokrovsky-Talapov type^ 
occurs and the ground state becomes a complete ordered 
ferromagnetic stated. 

The progress in the experimental front is achived 
by introduction of high-field neutron scattering studies 
and synthesis of magnetic quasi-one dimensional sys- 
tems. In many cases experimental data deviate signif- 
icantly from the theoretical predictions based on the 
pure isotropic AF Heisenberg model in the uniform 
magnetic field4i5i£iL2i2iI£. These deviations are due to 
anisotropies, most notably the Dzyaloshinskii-Moriya 
(DM) anisotrop y 11 ' 12 . The Hamiltonian of this model 
is written as 

6 = y*[j-§ j .-3 j+1 +{-iyv.(g j x-g j+1 )-HSf\, (2) 



where ~8 is the DM vector and the direction of this vector 
will be chosen along the y axis, 13 = (0,13,0). In actual 
systems, the direction of the ~l5 vector is fixed by the mi- 
croscopic arrangement of atoms and orbitals. Since the 
DM interaction breaks the fundamental SU (2) symmetry 
of the isotropic Heisenberg interactions, it is at the ori- 
gin of many deviations from pure Heisenberg behavior. 
Such anisotropy induces qualitatively different effects. In 
particular, in the ID AF spin- ^ model with the DM inter- 
action a gap is opened in the energy spectrum and scales 
as A ~ (DH)'s in contrast with the pure Heisenberg case 
(Eq.lJTJ). Theoretically, using bosonization techniques 
Oshikawa and Affleck^ explained the observed scaling 
behavior of the energy gap. They have shown that in 
the presence of the staggered DM interaction along the 
chain, an applied uniform field i? also generates an effec- 
tive staggered magnetic field h cx The staggered 
magnetic field for H — produces an AF ordered (Neel 
order) ground state and induces a gap in the spectrum 
of the model it is scaled as . 

For the higher-dimensional cases there is a theoreti- 
cal expectatio n 13 i 14 ' 15 that the field dependence of one 
of the gaps should be A ~ . Fouet et al. also 

studied— the gap-induced by the staggered magnetic field 
at the saturation uniform field H c = 2 J. Using field the- 
oretical arguments and density matrix renormalization 
group (DMRG) method, they found that the gap scales 
as A(H C ) ~ D$ . In a very recent work, this scaling be- 
havior is clarified by using exact diagonalization Lanczos 
results 16 . Also, it was shown that in the case of the 2D 
frustrated dimer singlet spin systems, a magnetic field 
induces staggered magnetization^. 

It should be noted that, most of the studies have ex- 
cluded from consideration the quantum effects associated 
with the DM interactio n 13 ' 14 ! 18 ! 19 ! 20 . In a recent work2i, 
the effect of an external magnetic field on the 2D AF 
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FIG. 1: The staggered chiral correlation function at H — 
plotted as a function of the DM vector D. The results 
reported for different chain lengths N = 12, 16, 20 and J = 
1. In the inset, the value of the staggered chiral correlation 
function C v is shown versus D for the pure DM interaction 
(H = and J = 0) and chain length N = 20. 



Heisenberg model with DM interaction is studied. The 
dependence of the quantum corrections on the DM inter- 
action is investigated. It is shown that the effect of the 
external field on the gap can be predicted by investigat- 
ing the on-site magnetization of the model. On the other 
hand, the interplay of DM interactions and an external 
magnetic field in spin-^ dimers is studied 2 ^. It is shown 
that the staggered magnetization of an isolated dimer 
has a maximum close to one-half the polarization, with a 
large maximal value of ~ 0.35 in the limit of very small 
DM interaction. They have also investigated the effect 
of the inter-dimer coupling in the context of ladders with 
DMRG calculations. However the interply of the DM in- 
teraction on the ground state properties of the ID AF 
spin-i Heisenberg model is much less studied. Since the 
integrability of the model will be lost in the presence of 
the DM interaction, very intensive studies are needed. 

In this paper, we present our numerical results ob- 
tained on the low-energy states of the ID AF spin-^ 
Heisenberg model with the staggered DM interaction (D) 
in an external uniform magnetic field (H). We study 
the mutual effect of a uniform magnetic field and DM 
interaction on the ground state phase diagram of the 
model. In particular, we apply the modified Lanczos 
method to diagonalize numerically finite chains. Using 
the exact diagonalization results, we calculate the spin 
gap, the magnetization, the staggered magnetization, the 
staggered chirality and various spin-structure factors as a 
function of the uniform magnetic field (H) and DM inter- 
action (D). Based on the exact diagonalization results, 
we obtain the ground state magnetic phase diagram of 
the model showing the Luttinger liquid, the staggered 



chiral, the spin-flop, and the ferromagnetic phases. We 
denote by "ferromagnetic phase" the phase with the mag- 
netization parallel to the external magnetic field as only 
nonvanishing order parameter. 

The outline of the paper is as follows: In section II we 
present our numerical results of the exact diagonalization 
calculations on the ground state properties of the AF 
Heisenberg chain with the DM interaction. In section 
III we investigate the effect of a uniform magnetic field 
on the ground state properties of the model. Finally we 
conclude and summarize our results in section IV. 



II. IN THE CASE OF H = 

In this section we explain the behavior of the model 
in the absence of the uniform magnetic field (H = 0). 
Classically, the effect of a staggered DM interaction is 
interesting. The DM interaction makes it energetically 
favorable for the spins to stay in the plane perpendicular 
to the direction of ~D {x-z plane). Without the "J" term, 
the DM interaction makes the spins in different sublat- 
tices to be under a 6 = ^ angle to each other. Thus at 
H = 0, in the ground state of the model (2) spins lie in 
the x-z plane (easy plane). 

On the other hand, Nersesyan et.al. predicted^ that in 
the anisotropyc (easy-plane) AF spin-^ chain with suffi- 
ciently strong frustrating next-nearest-neighbor coupling, 
a phase with a broken parity appears, which is charac- 
terized by the nonzero value of the chirality 



(3) 



where a denotes (%, y, z) and the notation (...) represent 
the expectation value at the lowest energy state. How- 
ever, two different types of the chiral ordered phases, 
gapped and gapless were found^^ 5 . In order to ex- 
plore the nature of the spectrum and the ground state 
phase diagram of the model, we have used the modi- 
fied Lanczos metho d 26 i 27 to diagonalize numerically fi- 
nite (N = 12, 14, ...,24) chains. The energies of the few 
lowest eigenstates were obtained for chains with periodic 
boundary conditions. 

We have calculated numerically the staggered chiral 
order parameter i Ylji^YXj an d the staggered chiral 
correlation function defined as 



1 N 



3 ^J+n/- 



(4) 



Lanczos results lead to \ — for any value of the DM vec- 
tor .D, because in a finite system no symmetry breaking 
happens. In Fig.l we have plotted the staggered chiral 
correlation function along the "y" axis, C v , as a function 
of DM vector D for different chain lengths N = 12, 16, 20. 
As can clearly be seen, the staggered chiral correlation 
function C y , increases with increasing D, which shows 
that the DM interaction suppresses the quantum fluctu- 
ations in the x-z plane and induces a staggered chiral 
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FIG. 2: The uniform magnetization M x as a function of 
applied magnetic field H for N = 20 chain for different values 
of the anisotropy DM vector D = 0.0, 0.1, 0.5, 1.0. The value 
of J = 1 is considered. 

phase in the ground state phase diagram of the model. 
Introducing a DM interaction, the SU(2) rotational sym- 
metry breaks and a quantum phase transition happens in 
the ground state phase diagram of the model. It is impor- 
tant to note that due to the profound effect of quantum 
fluctuations the chirality does not saturate. We have also 
checked the excitation energies of the three lowest levels 
as a function of D. We have considered the excitation 
gap in the system as the difference between the first ex- 
cited state and the ground state. We have found that a 
gap opens in the presence of the DM vector D. In the 
inset of Fig.l, the staggered chiral correlation function 
C y has been plotted as a function of D for the pure DM 
interaction (H = and J = 0) and chain length N = 20. 
It shows that in the absence of Hciscnberg interaction 
(J = 0) and Zeeman term (H = 0), the staggered chiral 
order is governed by D > 0. 

Thus, in the absence of a DM interaction, the ground 
state of the system is in the gapless Luttinger liquid phase 
with a power-low decay of correlations. Adding a DM 
interaction to the isotropic Heisenberg model develops a 
gap. The ground state then has the long-range staggered 
chiral order in the y direction. 



III. IN THE CASE OF H ± 

In this section we study the effect of a uniform mag- 
netic field on the ground state phase diagram of the ID 
AF spin-i Heisenberg model with DM interaction. As 
we mentioned, in the absence of the uniform magnetic 
field (H = 0) and DM interaction (D = 0), the spec- 
trum is gapless. The ground state is in the Luttinger liq- 
uid phase. By applying a uniform magnetic field H, the 



SU(2) symmetry of the pure Heisenberg model reduces to 
a U(l) symmetry corresponding to a rotation around the 
magnetic field direction (x axis). The spectrum remains 
gapless until a critical saturation field H c = 2 J. As soon 
as a staggered DM interaction with a i5 vectorc not par- 
allel to the uniform magnetic field (H) is introduced, the 
rotational symmetry in spin space is completely lost. The 
only symmetry that remains is the mirror symmetry with 
respect to the x-y plane (the plane contaning the uniform 
magnetic field and the i5 vector). As a consequence, the 
staggered magnetization per site must lie in the z direc- 
tion (the direction perpendicular to the plane defined by 
the uniform magnetic field and the ~D vector), while the 
uniform magnetization per site is in the field direction. 
If the 5 vector is parallel to the uniform magnetic field, 
the U(l) rotational symmetry is still present. Thus, the 
staggered magnetization is identically zero. In following 
we show that the quantum phase transitions can be eas- 
ily observed from the numerical calculations of the small 
systems. 

The symmetry breaking considerations suggest that an 
insight into the nature of different phases can be obtained 
by studying the magnetization 

MQ = ^£<^>< ( 5 ) 

j 

and the staggered magnetization 

M« = l^(-L)^}, (6) 

3 

and the spin correlation functions. The static spin struc- 
ture factor at momentum q is defined as 

S™(q) = J2e iqn (S?S? +n ). (7) 

n 

It is known that the spin structure factors give us 
deeper insight into the characteristics of the ground 
state. In particular we study the ii-dependence and in- 
dependence of different spin stucture factors. To deter- 
mine the properties of this model in different sectors of 
the ground state phase diagram we have implemented the 
Lanczos algorithm of the finite chains N = 12, 14, ...,24 
to calculate the lowest energy state. We have computed 
the ground state for different values of the DM vector 
D = 0,0.1,0.5,1.0. In Fig. 2 we have plotted the mag- 
netization along the applied uniform field, M x versus 
H for chain length N = 20 and different values of the 
anisotropy D. It can be seen, that for the D = 0, due to 
the effect of the quantum fluctuations in finite sizes the 
magnetization remains zero for smal values of the uniform 
magnetic field H . For H > H c = 2 J the magnetiza- 
tion saturate. This is in agreement with results obtained 
within theoretical approaches. Due to the quantum fluc- 
tuations, in the presence of D there is no sharp transi- 
tion to the saturation value of the magnetization. When 
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FIG. 3: The staggered magnetization along z axis, M* t as 
a function of applied magnetic field H for N — 20 chain for 
different values of the anisotropy DM vector D = 0.1, 0.5, 1.0. 
The value of J = 1 is considered. 



D / 0, the magnetization develops as soon as the mag- 
netic field is swiched on, only reaching saturation asymp- 
totically in the limit of infinite field. We mensioned that 
at the H = the staggered DM interaction causes the 
spins stay in the x-z plane. In this case, the effect of the 
uniform magnetic field decreases the degeneracy of the 
ground state energy. 

In Fig. 3 we have also plotted the staggered magneti- 
zation along the z axis M* t , as a function of the uniform 
magnetic field H. The results reported for a chain length 
N = 20 and different DM vectors D = 0.1,0.5,1.0. It 
shows in complete agreement with the theoretical results 
of the effective Hamiltonian^ and symmetry breaking 
considerations, by applying a uniform field H, a profound 
Neel order in the z direction induces. Which shows that 
there is long range spin-flop order along z axis. The os- 
cilations of M* t at finite N for smal values of DM vector 
D, are the result of level crossing between ground state 
and excited states of the model. There is also a maximal 
value for the staggered magnetization per site, around 
~ 0.3. This value is of the order of the maximal value 
of the isolated dimer, and it depends relatively weakly 
on D. The inset of Fig. 3 shows the staggered magneti- 
zation along the n z n axis versus the uniform magnetic 
field H at DM vector D = 0.5 and different chain lengths 
N = 12, 16, 20, 24. It can be seen that the maximal value 
of the staggered magnetization is independent of the sys- 
tem size. On the other hand, we have also investigated, 
S zz (q — 7r) as a function of the uniform field H for differ- 
ent DM vectors. We have found that, there is a trend to- 
ward staggered magnetization along "z" axis for H > 0. 
Which confirms that the ground state of the model has 
the Neel long range order along the "z" axis, which is 
known as the spin-flop phase. 
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FIG. 4: The value of the staggered chiral correlation func- 
tion along y axis, C y at the fixed different DM vectors 
D = 0.1,0.5, 1.0 versus the uniform magnetic field H for (a) 
J — and (b) J ^ 0(J = 1). The results reported for a chain 
length N = 20. 



An additional insight into the nature of different phases 
can be investigated by studying the staggered chiral cor- 
relation function C y . In this case we study the magnetic 
field dependence of the function C v for several values of 
the DM vector. In Fig. 4a we have shown the field de- 
pendence of the staggered chiral correlation function C y 
for the pure DM interaction (J — 0) and chain length 
N = 20. As we mensioned before, in the absence of 
the Heisenberg interaction and Zeeman term, the ground 
state of the model has the long-range staggered chiral or- 
der along the y direction. It can be clearly seen, that the 
staggered chiral phase remains stable even in the presence 
of the uniform magnetic field less than some critical field. 
Here a quantum phase transition occurs and the ground 
state becomes a completely ordered ferromagnetic state. 
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phases. The gapped staggered chiral and spin-flop phases 
are realized only in the case of DM interaction (D > 0) . 
Each phase is characterized by its own type of the long- 
range order: the ferromagnetic order along the magnetic 
field axis in the FM phase; the Neel order along the " z" 
axis in the spin-flop phase; and the staggered chiral order 
along the "y" axis in the staggered chiral phase. 

In principle, the DM interaction breaks the funda- 
mental SU(2) symmetry of the pure isotropic Heisenberg 
interactions and also U(l) symmetry of the Heisenberg 
model in the presence of a uniform magnetic field. The 
only symmetry that remains is the mirror symmetry with 
respect to the plane containing the uniform magnetic field 
and the ~l3 vector (x-y plane). 



IV. CONCLUSIONS 



FIG. 5: Schematic picture of the ground state phase diagram 
of the ID AF spin- i Heisenberg model with DM interaction 
in an external uniform magnetic field (H). 

Our numerical results also show that the value of the 
critical field depends on the DM interaction. 

To obtain a complete picture of the ground state phase 
diagram of the model, we have also calculated the func- 
tion C v for the Heisenberg model with the DM interac- 
tion. In Fig. 4b the staggered chiral correlation function 
is plotted as a function of the uniform magnetic field H 
for different values of the DM vector D = 0.1,0.5,1.0. 
The results reported for a chain length TV = 20. As is 
seen from this figure, C y first increases versus the uni- 
form magnetic field H and after passing a maximum de- 
creases. Since the uniform magnetic field suppress the 
quantum fluctuations in the x direction, thus the stag- 
gered chirality increases increases up to the saturation. 
However, with more increasing the magnetic field H, the 
staggered chirality decreases from the saturation value 
in well agreement with Fig. 4b. Also, due to the quan- 
tum fluctuations, the staggered chirality dose not reach 
to zero for large DM vectors. 

Based on the symmetry analyses and numerical calcu- 
lations, we expect that the ground state phase diagram of 
the model (2) has been a form shown on Fig. 5. Here the 
ground state phase diagram of a ID AF spin-i Heisen- 
berg model with the staggered DM interaction in an ex- 
ternal uniform magnetic field, is presented on the H-D 
plane. The ground state phase diagram contains, besides 
the gapless Luttingcr liquid and gapped ferromagnetic 
(FM) phases, the gapped staggered chiral and spin- flop 



In this paper, we have investigated the ground state 
phase diagram of the ID AF spin-i Heisenberg model 
with the staggered Dzyaloshinkii-Moriya interaction in 
an external uniform magnetic field H. We have imple- 
mented the modified Lanczos method to diagonalize nu- 
merically finite chains. Using the exact diagonalization 
results, we have calculated the spin gap, the magnetiza- 
tion, the staggered magnetization, the staggered chirality 
and various spin-structure factors as a function of the uni- 
form magnetic field and DM interaction. In the absence 
of the uniform magnetic field H = 0, we have shown that 
the DM interaction induces the staggered chiral phase. 
We have also found that a gap opens in the spectrum of 
the model. Thus, we have concluded that the gapped 
staggered chiral ordered phase appears in the ground 
state phase diagram of the ID AF spin-i Heisenberg 
model for D > 0. We have identified that the application 
of a uniform magnetic field induces a spin-flop phase in 
the ground state phase diagram of the model. Also, the 
staggered chiral phase remains stable even in the pres- 
ence of the uniform magnetic field. Finally, we have 
shown that the ground state phase diagram consists of 
four Luttinger liquid, ferromagnetic (FM), staggered chi- 
ral and spin-flop phases. 
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